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Abstract 

We study coherent pion production in neutrino-nucleus scattering in the energy region relevant 
to neutrino oscillation experiments of current interest. Our approach is based on a combined 
use of the Sato-Lee model of electroweak pion production on a nucleon and the A-hole model 
of pion-nucleus reactions. Thus we develop a model which describes pion-nucleus scattering and 
electroweak coherent pion production in a unified manner. Numerical calculations are carried 
out for the case of the 12 C target. All the free parameters in our model are fixed by fitting to 
both total and elastic differential cross sections for it— 12 C scattering. Then we demonstrate the 
reliability of our approach by confronting our prediction for the coherent pion photo-productions 
with data. Finally, we calculate total and differential cross sections for neutrino-induced coherent 
pion production, and some of the results are (will be) compared with the recent (forthcoming) 
data from K2K, SciBooNE and MiniBooNE. We also study effect of the non-locality of the A- 
propagation in the nucleus, and compare the elementary amplitudes used in different microscopic 
calculations. 
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I. INTRODUCTION 



The detailed theoretical study of neutrino-nucleus reactions is of great current impor- 
tance due to the ever increasing precision of neutrino oscillation experiments (recently car- 
ried out, on-going and forthcoming). Since most of these experiments measure the neutrino 
flux through neutrino-nucleus scattering, reliable theoretical estimates of the relevant cross 
sections are prerequisite for the accurate interpretation of the data. Some of these experi- 
ments (T2K, MiniBooNE, etc.) use neutrinos in an energy range within which the dominant 
processes are the quasi-elastic nucleon knockout and the quasi-free single-pion production 
through the excitation of the A (1232) resonance. Meanwhile, coherent single-pion produc- 
tion in this energy region (albeit not a dominant process) is also of considerable interest, 
since it allows us to study, with no ambiguity concerning the final nuclear state, the details of 
the A-excitation mechanism and medium effects on the pion; the knowledge of these details 
is essential for predicting the dominant quasi-free pion production processes. In this paper 
we focus on the coherent single-pion production process. 

There have indeed been quite active experimental efforts to investigate neutrino-induced 



coherent single-pion production in the A-excitation region. K2K[lJ and SciBooNE[2j investi- 
gated charged-current (CC) coherent pion production, while MiniBooNE (3] studied neutral- 
current (NC) coherent pion production. Furthermore, results for the anti-neutrino-induced 
coherent pion-production processes are expected to become available soon from MiniBooNE 
for the NC process 0|, and from SciBooNE for the CC process It is to be remarked, 
however, that the recent experimental results offer a rather puzzling situation. The exper- 
iments at K2K J and SciBooNE^ report that the CC process is not observed, whereas 
the MiniBooNE experiment (^J concludes that the NC process is observed. Now, from the 
isospin factors, we expect an approximate relation acc ~ Although the muon mass 

can reduce the phase space for the CC process at low energies, we still expect that acc 
should be of a significant size compared with <Jnc, and hence the above experimental results 
seem quite puzzling. In this connection it is to be noted that the MiniBooNE's use of the 
Rein-Sehgal (RS) model ^ in analyzing the NC data has recently be questioned 0: for a 
critical review of the RS model, see Refs. 0, 0]. The CC data analyses in Refs.jlJ 0] did 
not rely on a particular theoretical model for coherent pion production itself but, in dealing 
with some other neutrino-nucleus reactions that entered into the analyses, certain models 
whose reliability was open to debate needed to be invoked. 

The theoretical treatment of coherent pion production can be categorized into two types: 
a PCAC-based model and a microscopic model. In the former approach, the hadronic matrix 
element for neutrino-induced pion production is related to the pion-nucleus (or pion-nucleon) 
scattering amplitude through the PCAC relation. Meanwhile, in the microscopic approach, 
the hadronic matrix element is calculated by summing the elementary amplitude for weak 
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pion-production off a single nucleoli embedded in a nuclear environment. 

A prominent example of the PCAC-based approach is the model due to Rein and Sehgal 
(RS model) . Because of its success in the high energy neutrino process ^ (E u > 2 GeV, 
where E v is the incident neutrino energy) and its simplicity, the RS model has been exten- 
sively used in analyzing data in neutrino-oscillation experiments. Several authors, however, 
have recently pointed out that the RS model does not give a reasonable description for rel- 
atively low-energy neutrino processes {E v < 2 GeV) 0, 0] , and that the use of the RS model 
may have led to the puzzling experimental situation currently facing us. There have been 
several proposals 0, [lo|, [ill, E3] to remedy some of the possible insufficiencies in the original 
RS model. 

Meanwhile, in order to build a quantitatively reliable microscopic approach, it is obviously 
of primary importance to start with a model that can describe with sufficient accuracy 
electroweak pion production off a free single nucleon. Furthermore, for pion production off a 
nuclear target, we need to consider medium effects such as the final-state interactions (FSI) 
between the outgoing pion and nucleus, etc. Recently there have been several microscopic 
calculation^ [13I . \\A. T15I]. the most elaborate one being that by Amaro et al. {?]]. These 



calculations differ in the way the elementary process {y^N — > fi + Nir) is modeled and/or 
in the way the medium effects are taken into account. For example, only the resonant A- 
excitation mechanism is considered in Refs. 0, Q], while the non- resonant mechanism is 
additionally considered in Refs. 0, [lH| • It was shown in Refs. 0, [3] that the inclusion of 
the non-resonant mechanism leads to a reduction of the cross section by a factor of ~ 2, 
even though both models are constructed in such a manner that the data for the elementary 
process are reproduced fairly well. 1 This result indicates the importance of modelling the 
elementary process with a sound and systematic approach which has been extensively tested 
by available data. 

The purpose of the present article is to develop an alternative microscopic model for 
coherent pion production. An important ingredient of our formalism is a reliable dynamical 



model for the elementary process, and for that we shall employ the Sato-Lee (SL) model 16, 



13]. The SL model was first developed as a systematic framework for studying the resonance 



properties by analyzing data on pion production in photon (electron)-nucleon scattering 
in the A-resonance region 16J, |l8|. The SL model treats the resonant and non-resonant 



mechanisms on the same footing, and is known to provide a reasonably accurate description 
of an extensive set of p ion production data. The SL model was further extended to the 
weak sector in Ref. [17], and was shown to be able to reproduce data for neutrino-induced 
pion production off a nucleon. As has been done in the previous microscopic calculations, 



1 Unfortunately, conclusive data for the elementary neutrino process are still lacking, which leads to 
theoretical uncertainty. 
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we also need to incorporate the nuclear medium effects. In the energy region of our interest, 
the A-hole approach has proved to be successful in describing various processes involving 
pion-nucleus dynamics. These situations motivate us to develop a model for coherent pion 
production by combining the SL model and the A-hole model, and this is what we attempt 
in this article. We shall limit ourselves here to a case where the target nucleus (and hence the 
final nucleus also) has spin 0, and employ a simplified A-hole model proposed in Ref. 19]. 
As for concrete numerical calculations, we concentrate on the 12 C target, which has been and 
will continue to be an important nuclear target in many of neutrino-oscillation experiments. 
To test the reliability of our approach, we first calculate observables for coherent photo-pion 
production on 12 C using the same theoretical framework and show that the calculated results 
agree well with data. We then proceed to calculate observables for coherent neutrino-pion 
production on 12 C and present numerical results that can be compared with the recent data 
from K2K and SciBooNE. We shall also present theoretical predictions for those quantities 
for which experimental data will soon become available. 

The fact that the previous microscopic calculations exhibit rather large model-dependence 
makes it particularly interesting to use the SL model, which has been highly successful in 
the single nucleon sector. The SL model provides a consistent set of amplitudes for pion 
production and pion-nucleon scattering on a single nucleon; all these amplitudes are obtained 
in a systematic manner from the same Lagrangian. In our approach this consistency can 
be further taken over to the description of the FSI between the final pion and the nucleus. 
Thus, based on the SL amplitudes, we can construct a pion-nucleus optical potential that 
is consistent with the transition operators for electroweak pion production off a nucleus. To 
the best of our knowledge, our approach is the first to provide a consistent framework for 
treating the medium effect on the pion and electroweak pion production on the same footing. 
This point is worth emphasizing because it is this consistency that enables us to predict cross 
sections for electroweak coherent pion production with no adjustable parameters, once we 
fix certain parameters (see below) relevant to medium effects by fitting to the pion-nucleus 
scattering data. 

Another point to be noted is that our model takes into account the non-local effect for 
in-medium A-propagation. For neutrino-induced coherent pion production, neither the RS- 
based nor previous microscopic models have included this effect. As pointed out in Ref. 20], 
the non-local effect could reduce the cross section by a factor of ~ 2 (~ 1.7) for E v = 0.5 (1) 
GeV. We consider it important to take due account of the possibly large non-local effect. 

Our calculation adopts the following procedure. We first construct a pion-nucleus optical 
potential, employing the SL ttN scattering (on- and off-shell) amplitudes as basic ingredi- 
ent. The medium modification of the A-propagation in a nucleus is considered with the use 
of the A-hole model 19j. All the free parameters in our model (spreading potential, phe- 



nomenological terms in the optical potential) are fixed by fitting to pion-nucleus scattering 
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data. After these parameters are determined, we are in a position to make prediction on 
the coherent pion production process. Before calculating the neutrino-induced process, we 
test the reliability of our model by comparing our predictions for the photo-induced process 
with data. After finding satisfactory results for the photo-process, we proceed to calculate 
neutrino-induced coherent pion production. 

The organization of this paper is as follows. Sec. [Til is dedicated to the explanation of 
our approach. We first introduce the elementary amplitudes of the SL model. We then 
give expressions for calculating the electroweak coherent pion production amplitudes in 
terms of the SL amplitudes and derive the cross section formulae. The expression for the 
constructed optical potential and its relation with the scattering amplitude are also given 
there. We present numerical results in Sec. 1111} and give a conclusion in Sec. IIVI Appendix IA1 
provides the definition of the multipole amplitudes, while Appendix [B] explains the Lorentz 
transformation used in our calculation. In Appendix O we give expressions for quantities 
that appear in the A-hole model. 



II. FORMULATION 

The kinematics of the reactions under consideration is as follows. We consider coherent 
pion production in neutrino(z/^)-nucleus(t) scattering: Vi{p v ) + t(p t ) — > £~(Pi) + n + {k) + t(p' t ) 
for the CC process, and vi(p u ) + t{p t ) — ► vi{j>'i) + n°(k) + t(p' t ) for the NC process; we 
also consider the antineutrino-counterparts. The four-momentum for each particle in the 
laboratory frame (LAB) is given in the parentheses. The four-momentum transfer from the 
leptons is denoted by = p^ — p'/ 1 . We choose a right-handed coordinate system in which 
the z-axis lies along the incident neutrino momentum p„, and the y-axis is taken along 
p v x p' e . In evaluating a nuclear matrix element, it is convenient to work in the pion-nucleus 
center-of-mass frame (ACM). The kinematical variables in ACM are denoted by qA, J^a, etc. 
We also work in the pion- nucleon CM frame (2CM), when calculating the elementary SL 
amplitudes. The kinematical variables in 2CM are denoted by q 2 , k 2 , etc. When working 
in ACM (2CM), we choose a coordinate system in which the z-axis lies along qA (^2) and 
the y-axis is along p U)A x p' e A (p„ )2 x p' e>2 ). 



A. The SL Model 



We express nuclear transition amplitudes for coherent pion production in terms of the 
elementary amplitudes derived from the SL model 17] ■ In this section, therefore, we intro- 
duce the SL amplitudes. The differential cross section in the LAB frame for pion production 
in the neutrino- nucleon CC reaction, ve(p v ) + N(pn) — > ^~(p'e) + n + {k) + N(p' N ), is given by 
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(cf. Eq. (10) of Ref. [17j) 

d 5 a = G 2 F cos 2 6 C f\k\ \k\ -k- (p v -p' L ) \^ \p' e \ \k\ 2 m 2 N L^W^ 
dE'^dtt^ 2 \u n E' N ) \p v \uo^E N E' N (2tt) 5 ' U 

where Gf = 1.16637 x 10~ 5 GeV~ 2 is the Fermi constant, and 8 C is the Cabbibo angle 
(cos 0c = 0.974). E' £ and u n are the energies of the final lepton and pion, respectively, m N is 
the nucleon mass, and En {E' n ) is the initial (final) nucleon energy. W^ v and L^ u represent 



the hadron and the lepton tensors, respectively, and their definitions are found in Ref. 17 ] 
[Eqs. (11), (12)]. The above cross section can be written as 

-i 



d 5 a G F cos 2 9 C f \ k \ \k\ — k ■ {p u — p[ 



dE' e dn' £ dn w 2 \u v E' N 





k\ 2 


(2tt) 5 


\Pu\ 



E'.2Pu,2 (2) 



x \ E Eir 2 ,(^ y -^)l 2 

S N S 1\T S B 



where F v and F are the transition amplitudes in which the hadronic vector and the axial- 
vector currents are respectively contracted with the leptonic current. The symbol (s' n ) 
is the z-component of the initial (final) nucleon spin, while s' e denotes the final lepton spin. 
The energies of the final lepton and the initial neutrino in 2CM are denoted by E' t 2 and 
p U) 2, respectively. F v and F A , including both hadronic and lepton currents, are calculated 
in 2CM, and then embedded in the cross section expression given in LAB. The factor F 2 l 
arises from the relevant Lorentz transformation (see Appendix B) : 



2L 



\ lUkEnE'n 



where u n ^, En,2 an d E' N2 are the energies of the pion, the incident nucleon and the final 
nucleon in 2CM. The spin structure of F v and F A can be parametrized as 

F v = —id ■ e±F^ — a ■ k 2 d ■ q 2 x e±F 2 — id ■ q 2 k 2 ■ €±F^ — id ■ k 2 k 2 ■ e±F^ 

-id ■ q 2 q 2 • eF^ - id ■ k 2 q 2 • eF% + id ■ k 2 e F^ + id ■ q 2 e F^ , (4) 

where e± = q 2 x (ex q 2 ) and 

F A = -id ■ k 2 a ■ e±F A - d ■ q 2 x e ± F A - id ■ k 2 d ■ q 2 k 2 ■ e ± Ff - ik 2 ■ e±F A 

-id ■ k 2 d ■ q 2 q 2 ■ eF A - iq 2 ■ eF 6 A + ie Q F A + id ■ k 2 a ■ q 2 e F A . (5) 



The lepton-current matrix element e M is given by e M = (£|^7 M (1 — r )^)'^ v \vt). We have 
introduced parametrization for F A simply via F A = d ■ k 2 F v . The amplitudes, F^ and 
F A , are expressed in terms of the multipole amplitudes E^,Mf± A ,S^{f~ and Lf ±1 which 
are functions of q 2 and W (the irN invariant mass) and computed in 2CM. Their explicit 
expressions are presented in Appendix lAl 
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In a coherent process on a spin-zero target under consideration, only the spin non-flip 
terms of the transition amplitudes contribute. We therefore can work with F v ^ defined by 

PV(A) = bftpVM} ; (6) 

where the trace is taken for nucleon spin space. Their explicit forms are 

F v = -k 2 ■ q 2 x e L F 2 v , (7) 

and 

F A = -ik 2 ■ e±Ff - ik 2 ■ q 2 k 2 ■ e ± F A - ik 2 ■ e ± Ff 

-ik 2 ■ q 2 q 2 ■ eF A - iq 2 ■ eF§ + ie F 7 A + ik 2 ■ q 2 e F A . (8) 

In particular, the resonant parts of the elementary amplitudes are given by 

F V R - F£ = (-2k 2 ■ q 2 x e ± M^ V - 2tk 2 • e ± E^ A (9) 
-Mk 2 ■ q 2 e S^ A + 4ik 2 ■ q 2 q 2 ■ eL { ^ A ) kf , 

where the suffix U R" stands for the resonant parts of the corresponding multipole amplitudes 
associated with the excitation of the A resonance. From the resonant amplitude we can 
factor out the A-propagator, D(W), as 

fZ-f£= n{ **>* ) , (10) 

R R D(W) 

and 

D(W) = W - m A - E A (W) , (11) 

where m A and Sa are the bare mass and self energy of the A-resonance, respectively. 

We next discuss the T-matrix element for ttN scattering, which serves as an input for 
constructing an optical potential for pion-nucleus scattering. A calculational procedure for 
the ttN T-matrix within the SL model can be found in Ref. 16|. A distorted wave obtained 



with this optical potential will be used to take account of the final-state interaction in 

coherent pion production. The T-matrix is decomposed into the resonant and non- 
resonant (t nr ) parts as 

A c ) _ Ac) Ac) /-,r,x 

where the superfix c specifies a channel; in our model the resonance amplitude exists only 
for the P33 channel. The on-shell component of the T-matrix given in Eq. f fT2l) is related to 
the phase shift by 

Ac) _ W e -1 
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where W is the invariant mass of the irN system, and u nj 2 = y^2 2 + rr 4 an d En,2 = 
\Jk2 2 + m 2 N are the on-shell energies of the pion and the nucleon in 2CM, respectively. 
The resonant amplitude is expressed as 

igw«, h; w) = _ WM^pW , (14) 

where F^na^) is the dressed 7riVA vertex, and D(W) is the A propagator introduced in 
Eq. (II ip . We note that the four-momenta, ki and fc^, are in general off-energy- shell. 

B. Coherent pion production in neutrino-nucleus scattering 

Similarly to Eq. (fl2|) for the irN scattering amplitude, the weak amplitudes F vlyA \ defined 
in Eqs. (EJ-fJHD, also have a resonant F^ A ^ and non-resonant F^J-^ parts. Accordingly, the 
transition amplitudes of coherent pion production on nuclei have the resonant and non- 
resonant parts. We now describe how these two components are calculated in our approach. 



1. transition matrix element: resonant part 

The main task in calculating the resonant part of coherent pion production on nuclei is 
to account for the medium effects on A propagation in the elementary resonant amplitudes 
F^ A \ Here we follow the procedure of the A- hole model of pion-nucleus reactions by 
modifying the A propagator in Eq. ffTTj) . Thus it is useful to first briefly explain how the A- 
hole model is formulated by considering the elastic pion-nucleus scattering; for a full account 
of the formulation see Refs. 



21 



22 



23 



24fl. 



The A- hole model is formulated within the projection operator formalism 211] . The nu- 
clear Fock space is divided into four spaces; Po, Pi, D and Q. The Po-space is spanned by 
the pion and the nuclear ground state, the Pi-space by the pion and one-particle one-hole 
states, the P-space by the one-A one-hole configurations, and Q = 1 — Po — Pi — D contains 
the reminder of the full space. A projected Hamiltonian is written as, e.g., Hp £> = P$HD. 
Starting with the Schrodinger equation in the full space {H\^) = E\fy)), we can apply the 



standard projection operator techniques [2 ll] to obtain an equation, defined only in the Po- 
space, to describe the pion-nucleus elastic scattering T-matrix. In the A-hole model, one 
further imposes the condition that the P-space is the doorway of the transitions between 
P = Po + Pi and Q spaces; namely Ppq = Pqp = 0. The pion-nucleus scattering amplitude 
due to the A excitation can then be written as 

Tp p (E) = Hp oD GAh(E)H DPo , (15) 

where the total energy defined in ACM (E + Am^) is given by 



E + Am N = q° A + Jq\ + {Am N f = Jk\ + + Jk\ + {Am N y , (16) 
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where A is the mass number. The A-hole propagator G^h in Eq. (|T5|) is defined by 



G A \ = D(E - H A ) - W el - E pauli - S spr . (17) 

Here D(E — H&) can be calculated from Eq. ffTTj) with being the Hamiltonian for the 
A-particle in the nuclear many-body system. The effects due to the Q-space are included in 
the so-called spreading potential, E spr . A microscopic calculation of the spreading potential 
is very complicated since it involves the calculation of pion absorption by two or more 
nucleons. It is therefore a common practice to determine £ spr phenomenologically by fitting 
to the pion-nucleus scattering data. Excitations to the Pi-space are included in the A self 
energy Ea(W) of D(E — H&) [see Eq. (II ip ] with a correction due to the Pauli blocking 
(Ep au ii). De-excitation to the Po-space is the rescattering in the elastic mode, and is denoted 
by W e i. in our actual calculation, we expand G^h in term of W e i, and the expansion series 
is resummed by solving the Lippmann-Schwinger equation. 

The calculations of the pion-nucleus scattering amplitude in Eq. (]15|) require a diagonal- 
ization of the A-hole propagator Gaji of Eq. (fTTj) . For the diagonalization, it is practically 
convenient to work with the oscillator basis for the A state, defined by the Hamiltonian 
and the nucleon hole state. This diagonalization is a difficult numerical task. Although an 
efficient method using the doorway state expansion has been developed |23lj. the diagonal- 
ization of G/\h is still difficult, particularly for heavier nuclei. In Ref. [19|, Karaoglu and 
Moniz (KM) proposed a simplified calculation with the A-hole model in which G&u is calcu- 
lated with a local density approximation rather than a diagonalization. In their simplified 



treatment, E pau ij is calculated by a nuclear matter calculation 25|, and their result is given 
in Appendix Their parametrization of the spreading potential E spr in terms of a central 
and a spin-orbit terms are also given in Appendix Each term of the spreading potential 
has a complex strength, which are determined by fitting to the pion-nucleus scattering data. 
KM applied their approach to 7r- 16 scattering, and found a good agreement between their 



calculation with data, and also with the full A-hole calculation |22l . [23[ except for the most 
central partial waves. Encouraged by this success, we follow this simplified version of the A- 
hole model to include the medium effects on the A propagation in defining the electroweak 
pion production matrix elements. 

Schematically, the resonant part of the transition matrix element, of weak coherent 
pion production on nuclei induced by the charged current can be obtained by replacing the 
initial -Pdp of Eq. (fT5|) by Hdp^ where Pq is the space spanned by the (axial-) vector current 
and the nucleus in the ground state. In terms of the single particle wave functions i>j{pn) 
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of the nucleons in the initial and final nuclear states, we thus have 2 

KA A _ f d3pN d3 PN j.*r^' ^ 2A N(k2, q 2 ){2nf5(p N + q A - p' N - Ha) , , 
JVI R — 2s To ^ To \1™3 rVz? ; 7T^\ v v Yj\PN 

y J (27r) a (2iry J D(E + m N - H A ) - S pauU - S spr 

V- f d 3 p A , T 2A N(k 2 ,q 2 ) . , . 

= L / 7o~\3^ Pat) n , F , — v ^Wat) , (18) 

Y ^ ( 2n ) D \ E + m N- H A ) - Spauii - E spr 

where Pa = Pn + Qa = Pjv + ^a] the index j denotes single particle quantum numbers 
including the isospin . The summation is taken over the occupied states of the nucleus. 
The factor T 2A is defined by 



2A 



uJt,, 2 E n ^ 2 E' n ^ 2 
\ ^tt,aEn,aE' n ^ a 



where u n , and E' N are the energies of the pion, the incoming nucleon and the outgo- 
ing nucleon, respectively, and the quantities in the numerator (denominator) refer to 2CM 
(ACM). This factor arises from the fact that F^ A ^ computed in 2CM are to be embed- 
ded in Ai^ evaluated in ACM. To evaluate the numerator in the integrand of Eq. f[TB"j) . we 
clearly need a prescription for relating variables in 2CM to those in ACM. Here we use the 



commonly used prescription 27|, |28] to fix the nucleon momenta with the lepton momentum 



transfer q A and outgoing pion momentum k A as 

A-l , k A A — 1 

Pn = - ~YA~( qA ~~ kA > ' Pn = + ~Ya~^ Qa ~ A ' ' ^ ' 

and write the vrA^ invariant mass as 



W = ^{E NA + q\f - (p N + q A f , (21) 

with En a = \Jp 2 n + m jv- Having specified all the relevant variables in ACM, we can derive 
the corresponding variables in 2CM via a Lorentz transformation to obtain N(k 2 ,q 2 ) of 
Eq. (I18p . For more details about this Lorentz transformation (including the discussion of a 
somewhat different treatment of an off-shell pion momentum), see Appendix [B] Note that, 
in treating the wave functions, iP{pn) and xjj(p' N ), and the A kinetic term in the denominator 
in the integrand of Eq. (fTBl , we do not use the prescription given in Eqs. fTSOj) and (I2"T1) ; 
thus the important recoil effects on A-propagation are not neglected in our calculations. 
We incorporate the recoil effect on the A self-energy in the first orde r ap proximation. 



This is done by linearizing the A-propagator with the following expansion 2J] : 



In Rcf. 



D(E + m N - H A ) ~ D(W) - -y(W)(H A - e° A ), (22) 



the authors carried out a calculation for photon-induced coherent pion production by diag- 



onalizing Ga/i- 
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E + m N = W + e° A , (23) 
7(W) = 0Z>(W)/dW, (24) 

H A = ^ + V A + V% + e N , (25) 
1/fiA = l/m A + 1/(A - l)m N , (26) 

where V A (V a ) is the A (Coulomb) potential in the nucleus, and is the hole energy. 
The A potential is taken to be the same as that for the nucleon; its explicit expression 
is given in Appendix O Equation ( 1231) defines e A . To carry out the integration over the 
A momentum p A in Eq. (fl8l) . we express the nucleon wave function ipj(p) in terms of its 
coordinate-space form 4>j{r). We note that with the prescription in Eqs. (1201) and (1211) . the 
numerator N(k 2 , q 2 ) of Eq. (fT8l) is independent of the variable p A and can be factorized out 
of the integration. With this factorization approximation and with the use of the linearized 
form in Eq. ( |22l) . the integration over p A leads to the following r-space expression: 



M A R 



^T 2A N(k 2 ,q 2 ) \ ^ f ^y^v e^-j-l e ^ (f) f (2?) 

27T7 J ■ 1 7* 7*| 



where 



K\ = ^ {W - m A - E A (W) + 1 {E-W + m N )- 1 [e N + V A + ] - S pauli - E spr }(28) 



Following the procedure described in Ref. [19J [see Eqs. (25)-(39) therein], and subsequently 



applying the Lorentz transformation from ACM to LAB, we obtain the following expression 
for the transition matrix element M, R in LAB 



16^1 + |A|tt /UAj D(W) 1 



M L R = v ■ ■ ^ v > T2L (29) 
3 7 

x (-2fc 2 • q 2 x elMi 3/ 2' V - 2ik 2 ■ e ± E ( ^ A - Aik 2 ■ q^S^^ + Aik 2 ■ q 2 q 2 ■ eXg 7 ^ 



x E (! + ^r) / s 2 dsR 2 dRj ( P R)jo(Ps)^- (l + ^ [e* + p N (R)ji(k F s) 



where p = |fe^ — qr^l, P = | A*a + <Za|/2, s = |r' — r|, i? = |t* +7*|/2, and A"a is obtained from 
i^A by replacing with its average value, e^; we choose = 16 MeV. The 2CM variables 
k 2 and q 2 are obtained from k A and using the Lorentz transformation as mentioned above. 
The variable A denotes the charge state of the outgoing pion, while = 1 (—1) for N = 
proton (neutron). The factor T 2 l is from the Lorentz transformation from 2CM to LAB and 
is defined by 



2L 



^7r,2-E'Jv,2-E'i\r j 2 
\ ^tt,lEn,lE' n ^ l 



(30) 
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In Eq. ( I29l) . ji(x) is the spherical Bessel function of order £, and ji(x) = -ji(x); hp is the 
Fermi momentum 

k 3 F (R) = ^Pn(R) ■ (31) 

The proton (neutron) matter density is denoted by p p (p n ), and is normalized to the total 
number of protons (neutrons) inside the target. For the proton matter form factor we use 
the empirical nuclear charge form factor 29| divided by the proton charge form factor 30]. 
The neutron matter density is assumed to be the same as the proton matter density 
The single nucleon Hamiltonian appearing in Eq. f[2"9"l) is given by 



V 2 V 2 

' (R 2 + s 2 /A) l l 2 , (32) 



H N = -^5- - + V 



2 772 7v 8rriN 
where V is the single particle potential [Eq. (1C8I) ]. 

To take account of the final pion-nucleus interactions, we convolute the matrix element 
M L R of Eq. ([29]) with the pion distorted wave which is expanded in partial waves: 

X* x (k' A ) = E xiak' A )Y*^(k A )Y^(k' A ) , (33) 

/ 77 771 -n- 

where k' A is the off-shell momentum. We note that the pion distorted wave also depends on 
the pion charge (A). More details on our calculations of the pion wave functions are given 
in Sec. ITTDl 

By performing the partial wave decomposition of M R (now defined by the off-shell pion 
momentum by setting k A — > k' A ) and using Eq. fl33l . the amplitude M 1 ^ with pion-nucleus 
FSI takes the following form: 

M L R = e A J2 [PilM (cos0y^ - ism^l) (34) 



^71 



+Pi{x A ) (sin^ 2 + ioos^JfeJ) - 2P^x A )I l L ^ + 2P^(x A )I l s ^ 

where x A = q A ■ k A , (fi A is the azimuthal angle of the pion, and e A is the lepton current 
matrix element in ACM. The associated Legendre function of degree l n and order (1) is 
denoted by Pi n (Pi J)- We have introduced the quantities 1^ defined by 



= - t 32 y 1+ ^ MnfXA J dk ' Ak 1 xll ^ A ) J dx' A klY\ L Y 2ALl - l X R ^(x' A ) (35) 
E (1 + ^) / s 2 dsR 2 dR 3o (pR)j (Ps) e -^p- {l + ^ [e N + H N }} p N (R)Uk F s) , 



N=p,n 

where x' A = q A -k' A , x' 2 = q2-k' 2 , and 

21 



2k+±P e (x' 2 )P l7I (x' A ) , (X = L,S) 



12 



and 



X 



R 



D(W) 



p,(3/2),A Ayf (3/2),V r (3/2),A Q (3/2),A 
-^i? 1+ ' Jw i? 1+ > -^R 1+ 5 °i? 1 



+ ' 



(37) 



for X = E,M, L, S. 

The Lorentz transformation factors coming from the electroweak amplitudes (T2al) and 
the wave function (T x ) in Eq. f[3"4"|) are respectively 



2AL 



aE'n,lEn,l 



Utc,aE'x aE n A 



(38) 



where u;^ is the pion energy in the intermediate state, E 1 ^ and E^ N are the nucleon energies in 
the initial and final states while E' N and E'^ are those in the intermediate states; in general, 
E' N and E'^ can be different. As before, the suffices {2, A, L} attached to the energies specify 
reference frames. It is noted that the multipole amplitudes (X^) depend on x' A because the 
irN invariant mass in the intermediate state depends on it [Eqs. (1201) and (|21|) ]. We also have 
introduced the Lorentz matrix defined by e v 2 = A^e A ; also depends on x' A ; the same 
Lorentz matrix relates qA {k' A ) to q 2 (k' 2 )■ A procedure for deriving the Lorentz matrix and 
the transformation factors in Eq. (I38p are explained in Appendix [B] 



2. transition matrix element: non-resonant part 

We assume that there is no medium effect on the non-resonant part, F^ r — F^ r , of the 
weak pion production amplitude on a nucleon in nuclei. Including the final pion-nucleus 
interactions and using the same factorization approximation based on the choice Eq. (T2"U|) 
of the nucleon momenta to evaluate F^ r — F^ r , the non-resonant coherent pion production 
matrix element -M^ r can be written as 

Mt= E /^xl(^)r^r 2AL ^(^- qA )(F n ^-F^), (39) 

N=p,n 

where F^ (F^) is the non-resonant part of F v (F A ) given in Eqs. © and ©. F^ A) 
depends on N and A [Eq. (1A17h ] . and the set (A, A) is collectively denoted by Q. The nuclear 
form factor F^{p) is given by 

F N (p) = J d 3 rp N (ry p r . (40) 
After the partial wave expansion of the pion distorted wave, we arrive at 

M L nr = 6^ [p*(x A ) (cos&jy - isin^jfci) (41) 



I 7T 

+Pl(x A ) (sin^J^ 2 +ioos^jjjj) - Pi,(x A )jjrZ + PiAxa)J 1 s° 
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where we have introduced Jy " defined by 



rtr, = -47TZ J dk' A k'l x * xln (k' A ) J dx' A A;r* AL r 2AL 



'Xii — 
X 



E&to) E X% [ r 2 drp N (r)Upr) , (42) 

£ N=p,n J 

for X = E, M, L, S. The multipole amplitudes are included in as 

X £ J = {l + l?Xt\ + + fX^ e _ , (43) 

for X = L, S, and 

Eft = (i + l)E£f i+ - lEth- , (44) 
M% = {I + 1)M^ + + £M^_ , (45) 

for X = E, M. The ( dependence of the multipole amplitudes is indicated explicitly. For 
example, E^f g+ is the non-resonant part of Ef + which has been introduced previously. The 
same rule applies to the other multipole amplitudes. 



3. Cross Section 

Having written the transition amplitude for the coherent process in terms of the SL 
multipole amplitudes, we can proceed to calculate the cross section for the CC process. 
First, we write the transition amplitudes in Eqs. (134)) and (I4ip as 

M L R = M\^ A 
M L nr = M L m ,A ■ 

In the Laboratory frame, the differential cross sections for ue(p v ) + t(p t ) — > £~(pe) + n + (k) + 
t(p' t ) is then given by 



d 5 a G 2 F cos 2 9 C ( \k\ \k\ — k ■ (p„ — p' e ) 



dEidtt'AQ^ 2 {lUtt El 



x 

s 



t 

H |2 





Pill 


k 


2 


(2tt) ; 


''\Pu\ 



K,aPu,a (46) 



where E' t [= yp$ + (AmN) 2 j is the total energy of the nucleus in the final state in LAB, 
and E' eA and p u<A are the energies of the final lepton and the initial neutrino in ACM. Note 
that the calculation of \{Mr^ + M^)e A \ 2 of Eq. flU can make use of the following 
property: 

L A = 7, 2^ e A e A — Pu,APe,A + Pv,APe,A 9 Pv,A Pt,A Pv,ApPl,A° ' v 4 ' > 
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where g^ u is the geometric tensor and e^ upa is the antisymmetric tensor with e 0123 = 1. The 
plus (minus) sign in the last term is for the (anti-)neutrino process. 

To obtain the cross section formula for the neutrino NC process, v + t — > v + 7r° + t, we 
make the following changes in Eq. (|46|) : Remove the Cabbibo angle. Set the lepton mass 
equal to zero. Set the pion charge index A (and () to zero in I^" and J l ^ v (X = E, M, L, S) 
in Eqs. (1351) and (JUJ) . (Note that the pion wave function (xxi n ) a lso contains A-dependence.) 
Finally, multiply the multipole amplitudes M^! 2 ' 1 V with (1 — 2 sin 2 $w), where 8\y is the 
Weinberg angle (sin 9 W = 0.23), and multiply M$' v with (-2sin 2 w )- 

For the anti-neutrino CC process, the result for the neutrino CC process is modified as 
follows. Set the pion charge index A (and Q to —1 in I 1 ^ and Jjcu- Replace the lepton 
current by the one for the anti-neutrino process, which amounts to adopting the negative 
sign in the leptonic tensor, Eq. (I47p . What modifications are needed for getting the cross 
section for the anti-neutrino NC process is now obvious. 



C. Coherent Pion Photo-Production 



With the same derivation given above, we can also get an expression for the differential 
cross section of the coherent n° photo-production process 7(g) + t(p t ) — > ir°(k) +t(p' t ) in the 
LAB frame: 

where a is the fine structure constant, and \ J2 e stands for averaging over the photon polar- 
ization. The transition amplitudes M? R and M.^ r for the photo-process are obtained from 
Eqs. ( 134)) and (14T1) by retaining only the vector current, setting = 0, and regarding e\ as 
the polarization vector of the incident photon. Finally, the pion charge index (A) is set to 
zero in I 1 ^ and J 1 ^ [Eqs. (J35J) and 



D. Optical Potential for Pion-Nucleus Scattering 



We calculate the pion-nucleus scattering using the computer code, PIPIT (3l| by ap- 
propriately modifying the optical potential there to accommodate the dynamical features 
of the A-hole model and the SL model. In the original PIPIT, the optical potential (U), 
which is derived within the multiple scattering formalism by Kerman, McManus and Thaler 
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(KMT) [32J, is given by 



3 



A — 1 

fc A ) = — — {p p (qr)t 7rp (fc^, k A ; k° A ) + p n (q% n (k' A , k A ; k° A )} , (49) 

where fc^ (k' A ) is the incoming (outgoing) pion momentum in ACM, and k° A the magnitude 
of the on-shell momentum. The quantities p p (q) (p n (q)) is the form factor of the proton 
(neutron) matter distribution for q = k A —k' A , and t np (t nn ) is the pion-proton (pion-neutron) 
scattering T- matrix whose normalization has been defined in Eq. ( fl3l) . It is to be noted that 
this original optical potential does not take account of A-propagation in nuclei. In Ref. 19], 
KM separated t np (t^n) into the resonant and non-resonant parts, took the non-resonant and 
the Coulomb parts of the optical potential from the PIPIT code, and combined it with the 
resonant part derived from a simplified A-hole model. A phenomenological s-wave potential 
which is proportional to the square of the nuclear density (p t = p v + p n ) was also included 
to account for the pion absorption by two nucleons through non-A mechanisms. Thus the 
KM optical potential is given by 

U(k' A , k A ) = U nr + U R + U ph (p 2 t ) , (50) 

where U nr , Ur and U p h are the non- resonant, resonant and phenomenological parts, respec- 
tively. 

In constructing our optical potential, we follow the same separation as in Eq. fl50l . The 
non-resonant part of the optical potential is obtained from the PIPIT code by replacing the 
non-resonant T-matrices in the original code with those derived from the SL model. It is 
worth emphasizing that the SL model provides both on-shell and off-shell T-matrix elements. 
Another difference from the original PIPIT code is that we use a different prescription for 
the Lorentz transformation from ACM to 2CM, as explained in Appendix [B] 

Regarding the resonant part, we use the resonant part of nN T-matrix from the SL model, 
basically following the procedure used in Ref. [3] (apart from a more elaborate treatment 
of kinematics (Lorentz transformation, etc.)). First, we expand the optical potential into 
partial waves as 

U(k' A ,k A ) = - £ V l *(k' A , k A )Y£ mv (ti A )Y l7im% (k A ) • (51) 

7T , 



The resonant part of the potential is (cf. Eq. (39) of Ref. 191 ]) 
A - 1 8ir 2 p A 



x 



Vj?(V A , k A ) = — ^ J dx A T A2 - f - 1 x 2 P ln (x A )F 7TNA (k' 2 )F 7TNA (k 2 ) (52) 

E (1 + ^r) / s 2 dsR 2 dRj (pR)j (Ps)^ {l + ^ [e N + H N }} p N (R)Uk F s) , 



PIPIT also includes a finite-range Coulomb interaction, and corrections from the truncated part of the 



Coulomb interaction are taken into account using the Vincent-Phatak method 33 1 



1(3 



where k 2 (k' 2 ) is the incoming (outgoing) pion momentum in 2CM, and x A = kA ' k' A , 
X2 = k 2 ■ k' 2 , P = \k A + k' A \/2, p = \k A — k' A \. The dressed 7riVA coupling (-FVjva) has been 
introduced in Eq. f[T^|) . The Lorentz transformation of the T-matrix from 2CM to ACM 
gives rise to the factor Ya2 defined by 



A2 



^,2 u'k,2 E N,2 E N,2 



\ ^tt,A ^' n ,A En,aE' n ^ a 



with u& = JkP 2 + rnl, JJU = JkT + ml, E® = Jk^ 2 + m% and E$ A 



Pn,a + m N- The values of k$ and p$ A are fixed according to the prescription explained 
in Appendix [HI The other quantities have already been introduced in Sec. Ill B II 

Finally, we discuss the phenomenological term, U p h- We assume that in coordinate space 
U p h can be parametrized as 




u * r)=B m) ' (54) 

where B is the partial wave dependent strength of the potential. The corresponding partial 
wave potential in momentum space is given by 

Vf h (k' A , k A ) = ^4vr 3 ^ £ dx A P^x A ) J drr 2 Jo (pr) . (55) 

In the present calculation we include V p ° h and V p \ and treat their strengths B and B\ as 
adjustable parameters. Thus our model contains as free parameters B and B\ (complex 
numbers) in addition to the couplings in the spreading potential. 

Given the optical potential, we solve the Lippmann-Schwinger equation 

Ti(k' A , k A ; k A ) = VUk> A , k A ; k A ) + 2 -[ ^^j^^ ■ (56) 

ix J Un(k A ) + E t (k A ) — io n {k A ) — E t (k A ) + xe 

The solution to this equation will be used in two contexts. First, we use it to calculate 
pion-nucleus elastic and total scattering cross sections, and compare them with data to find 
the optimal values of the free parameters in our model. The solution to Eq. ( 1561) is also 
used to compute the pion distorted wave function that features in the matrix elements in 
Eqs. (1351) and (|42j) . For the former purpose, we obtain the full T-matrix of pion-nucleus 
scattering from T' in Eq. (1561) using the relation 

r . (57) 



A - I 



For charged-pion scattering, corrections for the finite range Coulomb potential are incor 



porated with the use of the Vincent- Phatak method 33J. The procedure for calculating 
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scattering observables from T is detailed in Ref. 31[- For the latter purpose, we calculate 
the pion distorted wave Xi (^a) associated with T' using the relation 

* ( u ) - 5 StA ~ k °A) | Tj(k%k A] k° A ) ( , R] 

XlA A) ~ k\ ^(k A ) + E t (k° A ) - uj n (k A ) - E t (k A ) + le ' 1 ) 

where, for notational simplicity, dependence on the pion charge (A) is suppressed. Following 
the KMT formalism 32|, we use xl (^a) i n evaluating the matrix elements in Eqs. (1351) and 



2J). This wave function is related to the full wave function by 

Xr^-^Y + ^XI- (59) 



For charged-pion scattering, xf * does not have the correct normalization, because the 
Coulomb potential has been cut off at a finite distance; this entails the necessity of multi- 
plying xf * with a normalization factor (call it k). We note that it is Xi rather than xlf 1 ^* 
that enters into our calculation, and we choose to use the same normalization factor k for 
Xi n as for Xj ■ Thus, in evaluating the matrix elements in Eqs. (13"5"j) and (l4"2"j) . we use kx^ 
instead of x* ■ m the A resonance region of our interest, it turns out that \k — 1| < 0.01. 
(For neutral pion scattering, k = 1.) 



III. NUMERICAL RESULTS 
A. Pion-Nucleus Scattering 

As explained in the previous sections, our model contains four complex free parameters. 
Two of them are the central (Vc) and LS (Vls) parts of the spreading potential [see Eq. (ICSp ]. 
and the other two are the strengths, B Q and Bi, of the s-wave and p-wave phenomenological 
terms in the optical potential [see Eq. ([5*5]) ]. These free parameters are optimized to fit the 
pion-nucleus scattering data. Since our aim here is to calculate coherent pion production 
off 12 C, we should use the it — 12 C scattering data to fix these parameters. Adjusting them 
to reproduce the total cross sections and the elastic differential cross sections for it — 12 C 
scattering, we obtain: 

V c = 48.0 - 34.5z MeV , V LS = -3.0 - 2.0i MeV 

B = 5.1 + 5.2i MeV , B 1 = 2.8 — 5.7i MeV . (60) 

We note that our calculations include the pion-nucleus partial waves up to l n < 9 [Eq. 056p ] . 
and s- and p-waves (and all possible spin-isospin states) for the elementary ttN scattering. 4 

4 Hereafter, we include the same set of partial waves (1^) in the amplitudes for both pion-nucleus scattering 
and pion production off a nucleus. For the non-resonant elementary pion production amplitudes, we 
include the partial waves up to £ < 4 in Eq. ([4*^|) . 
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Figures. Q] and [2] illustrate the quality of fit to the n — 12 C scattering data achieved in 
our model (with our optical potential). In Fig. [T], the total cross sections for 7r~- 12 C 
scattering are shown as a function of the pion kinetic energy T n in the laboratory frame. 
The results of our full calculation are given by the solid curve and, for comparison, the 
results obtained without the spreading potential are also shown in the dashed curve. We 
observe a large reduction in the total cross section as we go from the dashed to solid lines, 
which is mainly caused by the strong pion absorption simulated by the spreading potential. 
In connection with fitting to the pion-nucleus scattering data, it is worthwhile to make the 
following comment. In the calculation of coherent pion production, the final-state interaction 
is nothing but elastic scattering between the pion and nucleus. One might therefore think 
that a phenomenological adjustment of the pion-nucleus optical potential to fit the elastic 
pion-nucleus scattering data will be good enough. However, in our consistent model building, 
the spreading potential enters not only into the optical potential but also into the pion 
production operators, and hence it is important to control its strength using the total cross 
section data. The fact that the spreading potential has a very large effect on the total cross 
sections makes this point particularly important. 

Our results for the differential cross sections are shown in Fig. [2j In addition to our full 
calculation shown in the solid curve, we also show in the dashed curve the results obtained 
without the phenomenological term U p h [see Eq. (1541) ]. We see that this phenomenological 
p 2 term, which simulates absorption of s-wave and p-wave pions by two-nucleons within 
our model, is not large in the considered T v > 40 MeV region for tt— 12 C elastic scattering. 
However it is known that U p h can play an important role for many observables in low-energy 
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FIG. 2: vr- 12 C elastic differential cross sections. The solid curve is obtained with our full calculation 
while the dashed curve is obtained without the phenomenological terms in Eq. (|54p . The data are 
from Ref. for (a), Ref. 3 for ( b ) and Ref - H for (c)-(f). 



pion-nucleus scattering. As an example to shed light on this point, we have calculated 
7r — 16 elastic scattering at T n = 50 MeV using the same model (only the nuclear density 
is different). We have found that, in reproducing the data satisfactorily in our approach, 
the p 2 term plays an important role, its size being almost as large as that found in Fig. 4(a) 
of Ref. 



19( . Overall, the results of our full calculation satisfactorily reproduce the data for 



both the total and elastic cross sections. 



B. Coherent Pion Photo-Production 



We are now in a position to perform a parameter-free calculation of the cross sections 
for coherent pion production. The photo-process, for which extensive data are available, 
provides a good testing ground for checking the reliability of our approach. We compare in 
Fig.[3]our numerical resu 
the existing data 
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FIG. 3: Differential cross sections for 7+ 12 C 9 . s . — > 7r° + 12 C 9 . s . for different incident photon energies 
(indicated in each panel). The solid lines represent the results of the full calculation. The dashed 
lines are obtained without the FSI and without the medium effects on the A-propagation, while 
the dotted lines are obtained with the medium effects on the A included. The dash-dotted curves 
correspond to a case in which the pion production operator includes only the A mechanism. For 
more detailed explanations for the different cases, see the text. The data are from Ref. [37[ for (a) 
and from Ref. 



3J for (b)-(d). 



medium effects on A-propagation. 5 With the medium effects on the A included, the short- 
dash lines are obtained, and the results of our full calculation are given by the solid lines. 
Figure [3] indicates that the medium effects are quite sizable, and they play an important role 
in bringing the calculated differential cross sections in agreement with the data. Particularly 
noteworthy is the drastic reduction of the cross section in the A region [Fig. [3] (c)], a feature 
that reflects the fact that a significant part of the medium effects simulate pion absorption. 
The good general agreement seen in Fig. [3] indicates the basic soundness of the method we 

5 The "medium effects on the A" here refer to the combined effects of the Pauli blocking of A-decay 
(Spauii)j the spreading potential (S spr ), and the terms in the square bracket in Eq. (f^5|) . 
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have used in determining the spreading potential. 

It is true that, for higher incident energies, in the large angle region beyond the peak 
position, there are noticeable discrepancies between the results of our full calculation and 
the data. However, as noted in Ref . 38j] , the data in this region are likely to be substantially 
contaminated by incoherent processes in which the final nucleus is in its low-lying excited 
states. The effects of this type of contamination are expected to grow for higher incident 
photon energies and for larger momentum transfers (the large angle region) because of 
increased nuclear excitations. We therefore take the viewpoint that the discrepancy found 
in Figs. [3] (b)-(d) does not necessarily signal a failure of our model, and that our model 
describes coherent pion photo-production reasonably well. 

Figure [3] also shows (in the dash-dotted lines) the results corresponding to a case in which 
the pion production operator includes only the A mechanism (the non-resonant mechanism 
turned off); 6 the distorted pion wave function incorporating FSI is the same as that used 
for the full calculation. These results serve to demonstrate the importance of the non- 
resonant mechanism. Fig. [3] (a) indicates that, near threshold, the contributions from the 
resonant and non-resonant mechanisms are comparable, a feature that is not surprising away 
from the resonance peak. A remarkable feature is that even near the resonance energy [see 
Fig. [3] (c)] the contribution from the non-resonant mechanism is quite significant. This is 
partly because the resonant contribution is considerably suppressed by pion absorption (the 
spreading potential) and the non-local effect of A propagation (the A kinetic term). 7 

To summarize this section, the results for the coherent photo-pion production process 
establish to a satisfactory degree the reliability of our present approach (i.e., combined use 
of the SL model and the A-hole model) and motivate us to apply the same approach to 
neutrino-induced coherent pion production. 



C. Neutrino-Induced Coherent Pion Production 

We now present the numerical results of our calculations for neutrino-induced coherent 
pion production on the 12 C target. We consider the CC and NC processes induced by a 
neutrino or an anti-neutrino: 

Vp + 12 Cg. s . -> fT + 7T+ + 12 C 5 . S . 

v + l2 C g . s . -> u + 7i° + 12 C g . s . (61) 



6 In the SL model, the resonant amplitude itself contains the non-resonant mechanism. We refer to 
the purely non-resonant amplitudes as "non-resonant amplitudes" , and it is only these non-resonant 
amplitudes that we turn off here and later in Figs. and [TU1 

7 We come back to the non-local effect due to the A kinetic term later when we discuss the neutrino-induced 
processes. 
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^ + 12 C 9 . S . -> /i + + 7T- + 12 C 9 . S . 
* + l2 C,., - P + 7T° + l2 (',,,. 

Figure Hlgives the total cross sections for these processes as functions of the incident neutrino 
(anti- neutrino) energy in the laboratory system, E v . It is seen that, for higher incident 
energies, the ratio (JccI^nc approaches 2, a value expected from the isospin factor. For 
lower incident energies {E v < 500 MeV), however, a^c is larger than acc, reflecting the 
fact that the phase space for the CC process is reduced significantly by the muon mass. 
It is well known that interference between the vector and axial-vector currents can lead 
to different cross sections for the neutrino and anti-neutrino processes. However, since the 
coherent process is dominated by the contribution of the axial current (see Fig. [9]), the role 
of the interference term is diminished drastically. This explains why in Fig. H] the cross 
sections for the neutrino and anti-neutrino processes are almost the same. 

To compare our results with data, we need to evaluate the total cross sections averaged 
over the neutrino fluxes that pertain to the relevant experiments. We choose to use the 
fluxes up to E u < 2 GeV and neglect the fluxes beyond that limit based on the following 
consideration. Since our model includes no resonances other than the A, it is expected to be 
reliable only for W < 1.4 GeV. The fact that even at E v = 1 GeV coherent pion production 
can involve contributions coming from the W> 1.4 GeV region is disquieting, but we can 
still expect that the A-excitation contribution is predominant for the total cross section for 
the coherent process. [This feature can be seen in, e.g., Fig. [5] to be discussed later.] For 
E v ~ 2 GeV, we do expect that A dominance gets significantly less pronounced but that A 
still gives the most important contribution. Meanwhile, the region E v > 1.5 GeV belongs 
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to the tail of the neutrino flux used in MiniBooNE. We therefore consider it reasonable to 
compare with data our theoretical cross section averaged over the neutrino flux up to E u = 



2 GeV. For the CC process, we use the flux reported in Ref. 39J and deduce 

a™ = 6.3 x 10" 40 cm 2 . (62) 

A K2K experiment [l| reports the upper limit 

(TK2K < 7.7 x l(T 40 cm 2 . (63) 

In fact, this upper limit corresponds to events satisfying the muon momentum cut, p^ > 450 
MeV and the cut on the momentum transfer squared, Ql cc < 0.1 GeV 2 ; Q 2 CC is calculated as 

Q 2 TCC = 2E™{E I1 - p, cos 6,) - ml , (64) 

where the reconstructed neutrino energy (E r u cc ) is calculated from the muon kinematics [the 
energy (E^) and the scattering angle (6^)] assuming the quasi-elastic kinematics: 

rec = 1 (m 2 - m 2 ) - (m ra - Vf + 2E,(m n - V) 

2 {m n -V)-E^ + p^cos9 fl ' 1 ' 

where m p , m n and m M are the masses of the proton, neutron and muon, respectively and 
the nuclear potential (V) is set to 27 MeV. Our result in Eq.(l62l) is also obtained with these 
cuts, and is consistent with the K2K data. We note that a recent report from SciBooNE 
gives a similar empirical upper limit. 



For the NC process, we use the flux reported by MiniBooNE in Ref. 40] and arrive at 

= 2.8 x 10- 40 cm 2 . (66) 

This is to be compared with 

^MiniBooNE = 7.7 ± 1.6 ± 3.6 x lO^cm 2 , (67) 



given in Ref. (41J. Our result is consistent with the empirical value within the large experi- 
mental errors, even though the theoretical value is rather visibly smaller than the empirical 
central value. It is to be noted however that Ref. is a preliminary report, and that, as 
discussed in great detail in Ref. 0], cr MiniBooNE may be overestimated due to the use of the 
RS model fsj in the analysis. 

We now proceed to present our results for differential observables. In view of the fact 
that the event rates (cross section times flux) in the K2K, MiniBooNE and SciBooNE 
experiments 0, [3^] have been reported to have a peak around E u ~ 1 GeV, we shall often 
use this energy as a representative in the following presentation. Meanwhile, since the 
neutrino flux in the planned T2K experiment is expected to have a peak around E v = 0.6 ~ 
0.7 GeV 42J, we shall also present results for lower neutrino energies when that seems useful. 
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GeV; p-x is the pion momentum in the labo- 
ratory frame. The use of the solid, dashed, 
dotted and dash-dotted lines follows the same 
convention as in Fig. O 



FIG. 6: Same as in Fig. [5] but for E v 
GeV. 
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FIG. 7: The Q 2 -spectrum for za, 
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FIG. 8: Same as in Fig. but for the NC pro- 
cess at E,, = 1 GeV. 



The pion momentum spectrum for CC neutrino-induced coherent pion production is 
shown in Fig. [5] (Fig. [6]) for E u = 1 GeV (0.5 GeV). The importance of the medium effects 
manifests itself here in the same manner as in the photo-process (Fig. [3]). In the A region, 
strong pion absorption is seen to reduce the cross sections significantly, and FSI shifts the 
peak position. The dash-dotted line corresponds to a case in which the pion production 
operator contains only the A mechanism (without non- resonant contributions), while the 
pion optical potential is kept unchanged. We note that, at E v = 1 GeV (0.5 GeV), the dash- 
dotted line corresponds to 82% (64%) of the solid line (the results of the full calculation). 
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We have seen in the photo-process that the non-resonant mechanism is more important 
for a smaller energy transfer. To what extent the neutrino case should share this feature 
is not obvious because the axial- vector current contributions dominate here (see Fig. [9]). 
However, we can see in Figs. [5] and [6] that, in the neutrino case as well, the differential cross 
sections with smaller pion momenta are more enhanced by the non-resonant mechanism, 
and that this feature is more prominent for a smaller value of E v . A similar tendency is 
seen for the NC process also. These results indicate that the non-resonant amplitudes in 
our model, which are dressed by the rescattering, play a significant role in coherent pion 
production; their role is particularly important for E v < 0.5 GeV. This characteristic feature 
of our model should be contrasted with the fact that (tree-level) non-resonant mechanisms 
play essentially no role in any of the previous microscopic calculations for neutrino-induced 
coherent pion production. A more detailed comparison of the elementary amplitudes used in 
our present calculation and the previous microscopic-model calculations will be given later 
in Sec. UTTEl 

We show in Fig. [7J (Fig. [8]) the (^-distribution for the CC (NC) process. Note that Q 2 
defined by Q 2 = —q 2 = —{p v — p' e ) 2 is different from Q 2 ec defined in Eq. (jMj) . Because of 
the nuclear form factor effect, the distribution rises sharply as Q 2 approaches 0; for the CC 
process, however, (^-distribution becomes zero at Q 2 = due to the finite muon mass. Here 
again we show the results corresponding to a case in which the pion production operator 
contains only the A effect (with non-resonant contributions turned off). The non-resonant 
mechanism is seen to change the spectrum shape significantly and lead to a sharper peak. 

It is informative to examine the individual contributions of the vector and axial-vector 
currents. We show in Fig. [9] these individual contributions to the neutrino CC process. We 
find strong dominance of the axial- vector current. The nuclear form factor causes the drastic 
suppression of non-forward pion production. This aspect combined with the fact that the 
transverse photon coupling of the vector current [Eq. ([7])] forbids forward pion production 
leads to strong suppression of the vector current contribution. By contrast, since the vertex 
structure of the axial- vector current favors forward pion production, the strong suppression 
mechanism at work for the vector current does not apply here. This is the reason why 
the axial-vector current dominates. This result may be used to argue that incoherent pion 
production processes in which a nucleus does not break up but transits to excited states, 
are much less important than coherent pion production in the neutrino-nucleus scattering. 
As seen in Fig. El the incoherent processes give considerable contributions to the total 
pion production in the photo-process, 8 a feature that may lead to the expectation that 
the incoherent processes are considerable in the neutrino process as well. However, the 

8 The contributions from the incoherent processes are larger than they appear in Fig. [3] because sin 9^ 
needs to be multiplied in integrating over V . 
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FIG. 10: The effect of the non-locality of the 
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1 GeV. 



12 Cg.s.. The ratio 1Z of the total cross sections 
obtained with and without taking account of 
the non-local effect. 



mechanism responsible for the axial-vector dominance in the neutrino process works for the 
photo process in such a manner that coherent photo-pion production is strongly suppressed. 
Also, the inelastic transition form factor has a peak at a non-zero momentum transfer. As 
a result, for the photo reaction, the contributions from the incoherent processes become 
comparable to those from the coherent process. Thus the importance of the incoherent 
processes relative to the coherent process can be very different between the photo and 
neutrino processes. Takaki et al. 43| used a similar argument to explain a significant (very 
small) contribution from the incoherent processes in the photo-pion production (pion-nucleus 
scattering), compared to the coherent process. This argument may serve as a justification 
for the assumption currently used in data analyses that the incoherent processes need not 
be taken into account explicitly . 

Finally, we examine the effect of the non-locality of A-propagation in nuclei; because we 
employed the local density approximation for evaluating the A Green function [Eq. (117]) ]. 



. Although, 
201 . that study 



this effect arises only from the A kinetic term in the A Hamiltonian [Eq. 
as mentioned in the introduction, this subject has been studied in Ref. 
only included the A-mechanism without considering FSI or the medium effects on the A. 
It is thus interesting to revisit this problem in the framework of our significantly extended 
treatment. In the local approximation, we neglect the kinetic term in the A-Hamiltonian 
[Eq. (1231) ]. which means that the A is considered to be so heavy that it does not propagate 
in nuclear medium. To facilitate our discussion, we introduce the ratio TZ(E U ) denned by 

K{E V ) = a{E v )/a locA {E v ) , (68) 

where cr{E v ) represents the total cross section for z/ M + 12 C 9 . s . — > n~ + n + + 12 C„. S . calculated 



g.s. 
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with the A-propagator including the A kinetic term, whereas (T\ oca \{E v ) is that obtained 
in the local approximation. Figure [TU] shows TZ(E U ) calculated for the various cases. The 
long-dash curve corresponds to the A-only case (without FSI or the medium effects on the 
A; see footnote[5]) and the solid line to the case that includes the non-resonant components, 
medium effects on the A and FSI. To make comparison with Ref. 20[, we first consider the 
long-dash line; TZ(E U ) in this case is found to be 0.55, 1.03 and 1.14 at E v = 0.5, 1.0 and 
1.5 GeV. Meanwhile, Ref. Q reports K{E V ) = < 0.5, 0.6 and < 1 at E v = 0.5, 1.0 and 
1.5 GeV. Although both calculations indicate that the non-local effects are important, our 
results are qualitatively different from those of Ref. 20j. This difference originates from 
different ways of treating the energy in the A-propagator. In Ref. 20j, the in- medium A- 
propagator is assumed to be the same as the free A-propagator, whereas our A-propagator 
[Ga/i, Eq. (]17p ] is a nuclear many-body operator 24J (with some of the medium effects 
switched off). To illustrate this point, we include in Fig. [TU] (dash-dotted line) the results 
obtained with the use of the free A-propagator. In this case, we find 1Z{E U ) = 0.4, 0.76 
and 0.88 at E v = 0.5, 1.0 and 1.5 GeV, which is fairly close to the results in Ref. 20|. The 
result shown by the solid line indicates that, after the sophistication of the calculation, the 
non-locality due to the kinetic term is still important over the entire range of E v under 
consideration. In the previous microscopic calculations for neutrino-induced coherent pion 
production, the non-locality has not been explicitly taken into account. However, this does 
not necessarily mean that the earlier results are off by an amount suggested by comparison 
of the curves in Fig. [TO], for it is possible that the non-locality effects are partly included 
with the use of the spreading potential fitted to observables. In view of the importance 
of the non-local effect, however, we consider it preferable to take it into account explicitly, 
rather than include it operationally in the A mass shift. 

An additional point of interest is that it was reported in Ref. 20j that the non-locality 
changes the shapes of the differential cross sections. We remark that our results (not shown 
here) agree with that finding. 



D. Comparison with SciBooNE and MiniBooNE data 

The SciBooNE collaboration has been pursuing a further analysis of the data on neutrino 
and anti-neutrino CC coherent pion production, and some preliminary results have appeared 
in 

Refs.0, Q- 

These results contain detailed information on the differential observables 
for the pion and muon, and it seems informative to present our theoretical results in a 
manner that allows ready comparison with these data. To this end, we need to take into 
account the muon momentum cut (p^ > 350 MeV) and the momentum transfer cut (Qf. ec < 
0.1 GeV 2 ) adopted in the SciBooNE experiment; Q 2 rec has been defined in Eq. (1641) . The 
theoretical results we present in the following take account of these cuts unless otherwise 
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FIG. 12: The Q 2 CC distribution for 2/„ + 
12 C 5 . S . /i" + ^+ + 12 C 5 . S . at = 1 GeV. 



stated. We will present the results at E v = 1 GeV around which the event rate has a 
peak. Although, for direct comparison, we need to convolute the observables with the (anti- 
neutrino flux used in the SciBooNE experiment, the flux has not been released yet. We 
therefore present our results at a representative value of E v — 1 GeV. In Fig. [HQ, we show 
the cos ^-distribution for the neutrino and anti-neutrino CC processes. In the recent data 
analysis by the SciBooNE collaboration, events are classified according to the pion emission 
angle (9 n ). Their preliminary results exhibit a rather clear excess yield for 9 n < 35°, which 
is thought to be ascribable to coherent pion production. In our model, 85% of the pions are 
emitted in 9^ < 35° for the neutrino CC process at E v — 1 GeV, a feature that is in fair 
agreement with the preliminary SciBooNE result . 

Next we show in Fig. [T5] (solid line) the Q 2 ec distribution for the neutrino reaction. 9 
Only the cut is applied here for an obvious reason. We can see that the contribution 
from above Q 2 CC = 0.1 GeV 2 (the value adopted for the Q 2 ec cut) constitutes only a small 
fraction of the entire contribution (3% for the solid curve). The decomposition of the total 
contribution (solid curve) into two parts according to whether 9 n is smaller or larger than 
35° is shown by the dashed curve (9 n < 35°) and the dotted curve (9 n > 35°). The pion 
and muon momentum distributions are shown in Figs. [13] and [TH The upper (lower) end of 
the pion (muon) momentum distribution is sharply cut off because of the muon momentum 
cut > 350 MeV). The muon scattering angle distribution is shown in Fig. [TSJ Figures 
[TT]4T5] clearly show the characteristics of coherent pion production, i.e., sharply forward 
scattering (emission) of the muon (pion) with small momentum transfers. Finally, we show 



As discussed earlier, the neutrino and anti-neutrino cross sections differ only slightly. 
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FIG. 14: The muon momentum distribution 
for Vp + 12 C g . s . -» /i~ + 7T+ + 12 C 9 . S . at = 1 
GeV. 
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H~ + 7T+ + l2 C„. s . at E v = 1 



in Fig. [16] the spectrum with respect to the coplanar angle difference, A0, which is defined 
by A0 = (/v — 7r, where 0„- is the pion azimuthal angle in the LAB frame. (See Fig. [17] for 
a graphical representation of A0.) Fig. [16] shows slight asymmetry in the A<j) distribution 
around A0 = 0. It is interesting to note that this asymmetry is generated mostly by the 
contribution from the non- resonant amplitudes. To demonstrate this point, we present in 
the same figure the results obtained with the non-resonant amplitudes turned off, (dash- 
dotted curve). We also remark that the asymmetry arises mostly from the kinematical 
region satisfying 9 n > 35° (see the dotted curve). A similar asymmetry also arises for the 
anti-neutrino process. 

The SciBooNE collaboration have recently presented their preliminary results correspond- 
ing to Figs. [TTriTBI for both of the neutrino and anti-neutrino CC coherent pion production 
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FIG. 16: The coplanar angle difference (A0) FIG. 17: Graphical definition for the coplanar 
distribution for f A1 + 12 C 5 . s . — > [i~ +it + + 12 C g , s , angle difference (A0). 
at E v = 1 GeV. The definition of A(p is given 
in Fig. [T71 

reactions 0, Q]. When the flux prediction for the SciBooNE experiment becomes available, 
we will be able to convolute the results of our calculation with the flux and make direct 
comparison with the data. 

Meanwhile, the MiniBooNE collaboration has been investigating the NC process in 
(anti-) neutrino- nucleus scattering, and some results for the neutrino process have been pub- 
lished jsj , and more results are expected to be released. Since the neutrino flux information 
for the MiniBooNE experiment is available 40J], we can give the theoretical values of rele- 
vant observables convoluted with the flux. At present, data are publicly available only for 
the //-distribution [// = E n (l — cos #„-)], and we compare our calculation for this quantity 
with the data. In the analysis of the MiniBooNE NC data, the //-distribution was used to 
distinguish coherent pion production from other processes contributing to the 7r°-production 
events. To be more specific, MiniBooNE used the "shape" of the //-distribution obtained 
from the RS model jf| with the momentum reweighting function applied. It has been found, 
however, that a microscopic calculation in Ref. p| gives an //-distribution appreciably differ- 
ent from that obtained in the RS model, and the authors of Ref. 0] have pointed out that 
the MiniBooNE might have substantially overestimated the NC events. Figure [18] shows 
the "average" //-distribution obtained by convoluting the //-distribution given by our present 
calculation with the MiniBooNE neutrino flux 140J. For comparison, the figure also shows 



the MiniBooNE Monte Carlo results (cf. Fig. 3b of Ref. [sj), arbitrarily rescaled to match 
the theoretical curve at i] = 0.005 GeV. We remark that the //-distribution we have obtained 
is fairly close to that given in Ref. because the non-resonant amplitudes do not change 
the shape of the //-distribution significantly. Therefore, we arrive at the same conclusion as 
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in Ref. [7| that it is possible that MiniBooNE substantially overestimated the NC events. 

To facilitate a comparison of our calculation with data that are expected to be become 
available soon from MiniBooNE, we present theoretical predictions for some more quantities 
that are likely to be relevant. Figure [19] shows the flux-convoluted ir° momentum distri- 
bution predicted by our calculation. As far as observables for the anti-neutrino process 
are concerned, the flux-convoluted ^-distribution resulting from our calculation is given in 
Fig. [201 and the flux-convoluted tt° momentum distribution obtained in our model is shown 
in Fig. 



E. Comparison of Microscopic Models 

As mentioned, there are mainly two different theoretical approaches to coherent pion 
production in neutrino-nucleus scattering; a PCAC-based model and a microscopic model. 
The relation between the RS model (a PCAC-based model) and a microscopic model has 
been discussed in great detail in Ref. [7| , and comparison of those two models, including some 
improvement of the RS model, has been made in Refs. 0,8- The authors of Refs. 0,3 have 
emphasized that it can be problematic to use the RS model for E u < 2 GeV. To shed some 
more light on this issue, we consider it useful to make comparison of different microscopic 
models. In particular, we focus here on comparison between our model and the model of 
Wo et ah 0, which is the most sophisticated among the existing microscopic models for 
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neutrino-induced coherent pion production. The other microscopic calculations in the 
literature lack one or more aspects that are obviously important, such as the distortion of 
the final pion and the non-resonant mechanism for the weak currents. 

Here, we particularly focus on the elementary amplitudes for pion production off the 
nucleon. Our a ppr oach employs the SL model while Amaro et al. 0] used a model de- 



veloped in Ref. 46] (to be referred to as HNV). Both SL and HNV include the resonant 
and non- resonant amplitudes. A point to be noted, however, is that, although both mod- 
els reproduce reasonably well the data for the + N — > + 7r + + iV reactions after an 
appropriate adjustment of the axial- NA coupling, the two models involve rather different 
reaction mechanisms. In the SL model, we derive a set of tree diagrams from a given La- 
grangian with the use of a unitary transformation, and then we embed these tree diagrams 
in the Lippmann-Schwinger equation, which is solved exactly to yield a non-perturbative 
pion production amplitude that satisfies 7r-iV two-body unitarity. In HNV, on the other 
hand, a set of tree diagrams are calculated from a chiral Lagrangian. Then the sum of the 
contributions of these tree diagrams is identified with the pion production amplitude. At the 
tree level, the SL and the HNV models have essentially the same non-resonant mechanisms; 
a contact vertex in HNV may be interpreted as the vector meson exchange mechanism in 
SL. However the role of the non-resonant amplitude appears differently in the two models. 
In the SL model, non-resonant amplitude contributes constructively (destructively) to the 



10 A rather extensive comparison of numerical results from various calculations for the neutrino-induced 
coherent pion production, including those of Amaro ct al. Q, recent PCAC-bascd models and 
ours, has been presented at Nulnt09 by Boyd et al. [45 1. 
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resonant amplitudes below (above) the resonance energy. For + p — > pr + tt + + p, the 
interference of the non-resonant amplitude with the resonant amplitude changes in the SL 
model the total cross sections by a factor of 1.5, 1.02, 0.96 at E v = 0.5, 1, 1.5 GeV n , while 
the interference in the HNV always enhances the total cross sections; e.g., enhancement of 
a factor of 1.1 at E v = 1.5 GeV. The difference of the non- resonant mechanism appears also 
in the coherent pion production on 12 C, where only the spin and isospin non-flip amplitude 
contributes. Whereas the non-resonant amplitude plays an important role in our model (as 
seen in Figs. |5]and[6]), it plays essentially no role in the HNV model. In the neutrino CC 
coherent pion production, the full (tree) non-resonant amplitude increases the total cross 
section by 36% (19%) at E v = 0.5 GeV and 18% (0.4%) at E v = 1 GeV in our model. 
Thus the non-resonant mechanism in the spin-isospin non-flip amplitude is enhanced by the 
rescattering process. In the SL model, the non-resonant and resonant 7riV dynamics in the 
A resonance region has been tested using the extensive data of (7, 7r) and (e, e'ir) reactions. 
Although the SL model, which provides a unified description of the electroweak pion pro- 
duction reactions, describes very well the available data of the (v, in) processes, the current 
data do not yet allow to test the details of the reaction mechanism. 

Furthermore, utilizing the consistency of (y,£7r), (e, e'ir) and (it, it) reactions in the SL 
model, we have developed a model which treats photo- and neutrino-induced coherent pion 
production processes in a unified manner. Thus we were able to calibrate the reliability of 
our approach with data for the photo-processes, which is an aspect specific to our approach. 

IV. CONCLUSIONS 

We have developed a microscopic dynamical model for describing neutrino-induced co- 
herent pion production on nuclei. Because experimental data for neutrino (both elementary 
and nuclear) processes are rather limited, it is not straightforward to assess the reliabil- 
ity of theoretical calculations. A reasonable strategy to take seems to develop a model 
which describes strong and electroweak processes in a unified way, and then to test the 
model extensively by comparing with a large collection of data for the strong-interaction 
and photo-induced processes and with limited available data for weak processes. We have 
carried out this program here for the case of the neutrino-induced coherent pion production 
process. By virtue of the mentioned strategy, our model is probably the most extensively 
tested one among the existing models for this process. To achieve the stated goal, we need 
a theoretical framework that provides a unified description for the elementary (jr, ir'), (7, n) 
and (u, £ir) processes on a single nucleon. We have adopted the SL model, which is known 
to give satisfactory descriptions of these elementary amplitudes. We then have combined 

11 See footnote [SI 
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the SL model with the A-hole model to construct a theoretical framework that can describe 
in a unified way pion-nucleus scattering and electroweak coherent pion production. The 
unified nature of this approach allows us to fix free parameters in the model using the data 
for pion-nucleus scattering, which in turn enables us to make parameter-free predictions on 
electroweak coherent pion production off a nucleus. Another benefit of the present unified 
approach is that we can assess the reliability of our model by comparing the results for coher- 
ent pion photo-production with data. Our model is found to describe reasonably well both 
pion-nucleus scattering and coherent photo-processes, which establishes a basis for applying 
the same model to the neutrino-induced processes. 

Comparing our numerical results with the recent data on neutrino-induced coherent pion 
production, we have found that the result for the CC process is consistent with the upper 
limit from K2 K Q, and that the result for the NC process is somewhat smaller than the pre- 
liminary experimental value from MiniBooNE 4l| . However, as discussed in the literature, 
MiniBooNE's analysis may have overestimated the cross section due to the use of the RS 
model in their analysis. We have examined to what extent the various aspects of physics 
involved in our model individually affect the cross sections. We have shown that the medium 
effect on the A (the spreading potential effect in particular) and the FSI change the cross 
sections significantly. It is to be noted, however, that these rather drastic changes in the 
cross sections due to the medium effects are well under control because: (i) the spreading 
potential and the pion distorted wave function have been fitted to and tested by the empir- 
ical total and elastic cross sections for pion-nucleus scattering in and around the A region; 
(ii) the medium effects of a similar magnitude for the photo-process have been shown to 
bring our calculation into good agreement with the data. 

An interesting feature of our model is that the unitarized non-resonant amplitudes give a 
significant contribution to the cross sections. This is in sharp contrast with the results of the 
previous calculations; for instance, the calculations in Refs. 0,0], which considered a tree- 
level non-resonant mechanism, found almost no contribution from it. It is worth emphasizing 
that this noticeable difference should not be taken as a measure of uncontrollable model 
dependence because (as we confirmed) the difference arises largely from unitarization of the 
non-resonant amplitude, which clearly needs to be implemented. 

We have reexamined the non-local effect in A-propagation in nuclei. It was emphasized 
in Ref. 20[ that this non-local effect, despite its large size, was not considered explicitly in 



any of the existing models for neutrino-induced coherent pion production (whether based 
on a microscopic model or the RS model). The authors of Ref. [2^] made this remark 
based on their calculation that only included the A mechanism. Our present calculation, 
which additionally incorporates the spreading potential and FSI, also indicates that the non- 
locality gives a large effect. Thus, regardless of the level of sophistication in the treatment 
of medium effects, one should always include the non-locality effect explicitly. 
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Because it is expected that the SciBooNE and the MiniBooNE collaborations will report 
more detailed data on (anti) neutrino-induced coherent CC and NC pion productions, we 
have presented numerical results relevant to these experiments. 

Finally, we made a com par ison of the elementary amplitude (HNV46J) used by Amaro 



et al.[7| and ours (SL [ly, |17||) to clarify similarities and differences between them. The 



noteworthy points are: (i) At tree-level, both SL and HNV have essentially the same non- 
resonant mechanism; (ii) In the SL model, a unitary pion-production amplitude is obtained 
by solving the Lippmann-Schwinger equation in which the tree-diagrams are embedded, 
whereas, in the HNV model, the sum of the tree-diagrams are identified with the pion- 
production amplitude; (iii) The non-resonant amplitudes of SL and HNV work differently 
both for the elementary processes (e.g., -\- p —* /i~ + ir + + p), and for coherent pion 
production; (iv) In SL, the rescattering contribution contained in the non-resonant amplitude 
considerably enhances the cross section for coherent pion production. 
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APPENDIX A: MULTIPOLE AMPLITUDES 



The amplitudes , F^ in Eqs. (0J and ([5]) are expressed in terms of multipole amplitudes 



EYl A , M,T, SrT and L& as 



r V,A Q V,A 



J l± i lvl l± i u l± 



l± 



F? = T\ P U E l + PUEl + IPUK + + l)PLxM] 

= E^,^ + p'UEl - p;' +1 Mf + + piumI] 

I 

fJ = T\-Pi E l - p'i' E l + p"K - p" m I] 



F^Ttiil + VPULl-lPUL 
I 



V 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
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and 



K = £["(' + + lP(Ll\ , (A6) 

i 

= EhC + !)^+ + ^-1 > ( A7 ) 
i 

Fs = EK' + " ^-i^-] > ( A8 ) 

^ = + P'iFl + (I + 2)P[M* 1)P/M£] , (A9) 

F 2 A = EK' + + ^-i^-l , (A10) 

# = E^"< + P\ F\L + P'{Mt - P'{Ml\ , (All) 
i 

Ft = Et-^i< - - + K-iM£\ , (A12) 

^ = EI-C + + ^-1 - (A13) 

# = EK* + " lP ULf-\ , (A14) 

^ = EK' + VP»iS& - iPUst] , (Ais) 
i 

F S A = EH 1 + l ) P K + IPlStl (A16) 
i 

Pl{x) is the Legendre function and x = k ■ q; k and q are the pion momentum and the 
momentum transfer to the nucleon, respectively. 

The multipole amplitudes from isovector currents are further decomposed according to 
the total isospin (T) in the final ttN state as 

x£*= E *L T)V 'X-, (ait) 

T=l/2,3/2 

with X being E, M, L or S. We have introduced the projection operator Ajj defined by 

3/2 _ 2djj — itijk T k 

1/2 _ Qjj + iejjkTk 
3 

where the indexes i and j refer to the final pion isospin state and the component of the 
isovector current, respectively. For electromagnetic or NC processes, M^ V Ti, which is due 
to an isoscalar current, is also added to Eq. (1A17I) . 

In the main text we sometimes use the notation X]± A ^'^ , where ( collectively denotes the 
pion charge and the nucleon isospin state; is a matrix element (in isospin space) of 

Eq. (1A17I) . Since we are only concerned with coherent pion production, the specification of 
the pion charge determines % and j in Eq. (1A17j) . We can find the matrix element (in isospin 
space) of Eq. (IA17I) by specifying the nucleon isospin state. 



A 3/2 = toij Kijk'k (Alg) 
A l/2 = u i,j T ^ijk 1 k ^ (Alg) 
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APPENDIX B: LORENTZ TRANSFORMATION FROM ACM TO 2CM 



In coherent pion production in neutrino-nucleus scattering {yi + t — > i~ + tt + + t), the 
elementary process is W + (q A ) + N{p^) — > ir + (k A ) + N(p' N ), where the four-momenta in the 
pion-nucleus center-of-mass frame (ACM) are given in the parentheses. We suppose here 
that the pion momentum is on-shell. In a prescription we employ, the nucleon momenta are 
fixed as 

<1a A-l , k A A- 1 

Pn = --j - (Qa - k A ) , p N = -— + -^-{Qa - k A ) , (Bl) 

and the invariant mass (W) of the pion and nucleon is 



W = ^{p% + q A Y - (p N + q A f , (B2) 

where p% is the nucleon energy on the mass-shell. We note that W depends on x A {= k A ■ q A ) 
as well as \q A \ and \k A \. For convenience, we write W(|<7a|, I^aI; x a)- 

We perform the standard Lorentz transformation from ACM to the ttN CM frame (2CM). 
An arbitrary four-momentum in 2CM (p 2 ) is written with the corresponding four-momentum 
in ACM (p A ) as 

P2=PA~^P+ —^(PA ■ P)P , (B3) 



Pi 



P° Pa - Pa -P 



W 

with P = p N + q A . 

We now consider a case in which the pion momentum is off-shell (k' A ). We encounter this 
situation when we consider the final-state interaction in the coherent process. As before, 
the nucleon momenta are fixed using Eq. fIBip with k A replaced by k' A . However, we do not 
use the nucleon energy on the mass-shell. Instead, we take p° N so that 

W(\q A \, \k' A \,x' A ) = W(\q A \, \k A \,x A ) foix' A = x A , (B4) 

where W is obtained with Eq. (1B2I) . With the nucleon four-momentum (pn) obtained 
in this way, we can perform the Lorentz transformation as Eq. (IB3I) . This prescription 
greatly reduces the amount of labor involved in our numerical calculation, because the SL 
amplitudes need to be calculated at each value of W. With the variables obtained above, 
we can calculate T 2 al used in Eqs. (l3~5l) and (|42l) : 



2AL 



UJ 'n,2P'N,2PN,2 /r, c \ 

(B5) 



\ ^,aPn,lPn,l 



with u' A = Jk' A + ml. 
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Finally, we discuss the factor r x ,used in Eqs. ( 1351) and (|42|) . which originates from the 
pion wave function due to the Lorentz transformation. Among the final-state interactions, 
the simplest process is the scattering of the pion off a single nucleon 7r(k' A ) + N(p'^) — > 
+ N(p N ), where the variables in ACM are shown in the parentheses; only k A is on- 
shell. Similarly to Eq. (IBip . we fix the nucleon momenta as 



k' A A- I 



k a k / 



f k A 



A-l 



k a k f 



(B6) 



2A v A 1 ri " A 2A 

We assume here that the energies of all the nucleons are on the mass-shell. For the Lorentz 
transformation from ACM to LAB specified this way, we can calculate the Lorentz factor as 



U-k,aE'n,aEn,A ^ 



(B7) 



Although the actual final-state interaction includes multiple scattering processes, it is beyond 
our framework to calculate T x with multiple scattering taken into account. We therefore use 
r x calculated for the elementary process in Eqs. ( l35l) and ( l42l) . Actually, the Lorentz factor 
for the plane wave term in Eq. (|58|) is given by the the rightmost expression in Eq. (IB 71) . 
Because the approximate equality in Eq. (1B7I) is quite accurate for k' A = k A , we use the 
middle expression in Eq. ( 1B7I) to evaluate the matrix elements in Eqs. ( 1351) and ( l42l) . 

APPENDIX C: EXPRESSIONS FOR SOME COMPONENTS IN THE A PROP- 
AGATOR 



1. Pauli correction to the A self energy 



We follow Ref. 



25| to calculate the Pauli correction to the A self energy (Sp au ii). The 



irNA coupling is from the SL model. 



m N 



J Pauli 



w 

k F +(3 



k F ~l3 



26(k F - 0) I dqq 2 ^— 

Jo K z — q z + le 



(CI) 



q 2 + P 2 ~ k 2 F \ cu n (q)F^ e A (q)F nNA (q) 



dqq A 1 



i\k F -p\ \ 2q/3 J K 2 — q 2 + ie 

where 6(x) is the step function, kp is the Fermi momentum [Eq. (13TT) ]. W is the irN invariant 



mass [Eq. (T2TT) ]. co n (q) = \Jq 2 + m 2 , and 

m N r 



K 



W 



(W - m N 



,2 2 

— 



Furthermore, for electroweak pion production amplitude [Eq. (|35|) ]. 

m N 



w 



(Pn + qjC, 



(C2) 



(C3) 
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where Pn is fixed using Eq. (|20|) . and qA is the momentum transfer to a nucleus in ACM; 
for the optical potential [Eq. f)52p ]. is replaced with (the incoming pion momentum). 
We use the on-shell pion momentum to fix p^. The dressed ttNA vertex (-F^ata) is taken 



from Eq. (|T4|) . and the bare 7riVA vertex denoted by F^ A j s gj veri as |xg| 



7 bare / \ _ • fnNA 



E N (q)+m N ( Al NA 



\ 247r*E N (q)u„(q) \Al NA + q' 



(C4) 



2. A spreading potential 



We consider the following spreading potential consisting of the central and the LS parts: 

S spr = Vc^r + V LS f L s(r)2L A ■ E A , (C5) 
hs(r) = pr\~^ , (C6) 

with fi = 0.3 fm~ 2 . We have two complex coupling constants Vq and Vis which are fitted to 
pion-nucleus scattering data. The radial dependence of the LS spreading potential is taken 
from Ref. 47(. We implement the spreading potential [Eq. (105[) ] in the A-propagator after 



evaluating the doorway state expectation value of the LS term. Thus, the LS term provides 
an L-dependent shift of the resonance mass and width asj^] 

V L k T/ (0L\ptfLsk 2 - (pths)' i + (Pths)' \4>l) (nrA 

^LS - ~ bV LS J—, J- 2 , \i d \l \ ' y^') 

(<pL\ptk 2 ~ {pt)' Tr \<pL) 

with the plane wave pion function 0L(r) = ji{kr). 
3. A (nucleon) potential 



V A (r) = V(r) = (-55 MeV) . (C8) 



10 



4. A Coulomb potential 



(r > R e ) 



(r < R e ) 



(C9) 



f 2(Z-l)a 



Za 



Za 

T ' 



(Z- l)ar 2 , 3(Z- l)a 



Zar 2 i 3Za 

o , 

Zar 2 _ 3Za 
~2RT 2R e ' 



(tt + + n -> A+) 
(tt-+ P ^ A ) 
(7r~ + n -> A") 



In the above Z is the atomic number. The equivalent square well radius, denoted by R e , is 
related to the mean square radius ((r 2 )) of a nucleus by 



Rp 
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